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Enhancement of the trapping and optimization of the beam quahty are two key issues of Laser 
Wake Field Acceleration (LWFA). The influence of stochastic acceleration on the trapping of elec- 
trons is compared to the one of cold injection. It is shown that when considering a high intensity 
wave perturbed by a low intensity counter-propagating wave, in the non-linear blowout regime, 
the influence of the colliding pulses polarizations (either parallel linear or positive circular) on the 
beam quality seems weak when the electron density is below ^ 10~^ critical density. The effect 
of a homogenous constant magnetic field Bq, parallel to the direction of propagation of the pump 
pulse, is studied in the blowout regime. Transverse currents are generated at the rim of the bubble, 
which results in the amplification of the Bo field at the rear of the bubble. Without Bo field the 
beam periodically explodes and re-confines, this phenomenon is suppressed when Bo reaches some 
threshold, which is a function of the laser pulses parameters (intensity, waist, duration). Therefore 
the dynamics of the beam is modified, its maximum energy is slightly boosted and above all trans- 
verse emittance reduced. Moreover the low energy tail, observed in the non magnetized case, can 
be completely suppressed leading to very sharp mono-energetic beam when Bo is applied. If the 
available Bo field is limited then one has to fine-tune the spatio-temporal shape and intensity of the 
colliding pulse in order to get an acute control on the beam quality. 



INTRODUCTION 



In laser- Wakefield acceleration (LWFA) [H-Q, a laser 
creates a plasma wave wakefield with a phase velocity 
close to the speed of light (c). The acceleration gradi- 
ents in these wakefields can easily exceed 100 GeV/m, 
hence a cm-long plasma based accelerator can produce 
GeV-energy electron beams. A particle injected in such 
a wave gains energy from the longitudinal component of 
the electric field, as long as the pump pulse is not de- 
pleted and the dephasing length is not reached. These 
wakefields have ideal properties for accelerating electrons. 
The transverse focusing field increases linearly with the 
radial distance and the accelerating longitudinal field is 
independent of the radial coordinate [5, 6]. LWFA can be 
split into three different options . The first corresponds 
to a plasma density rie ^ lO^^cm"^, a pulse length (cr) 
matching half of a plasma period and a spot size (wq) 
roughly equals to the bubble radius, wq ^ cr ^ 
where ao is the normalized vector potential of the laser. 
This is the idea of the bubble regime 0, Q . For these 
conditions, a hundred-joule class laser would have an 
intensity of the order ~ 10^^ W cm~^. In this regime, 
the electrons are continuously injected, this results in 
tremendous beam loading and the loaded wake is noisy. 
In this paper we explore different techniques to improve 
the beam quality of LWFA when electrons are injected 
in the wake with a colliding pulse. Hence, the bub- 
ble regime is not appropriate. We rather select mod- 
erate laser intensity / < lO^^W.cm"^ and plasma den- 
sity Tie < lO^^cm"^ according to the guidelines proposed 



by Lu et al.[9] to achieve a more controlled and stable 
blowout of the electrons. Self-injection of electrons can 
occur when the pump pulse intensity is high but the ac- 
celerating structure is almost the same. In this frame 
beam loading effects clamp further injection leading to 
beams with a smaller energy spread. In order to limit 
the computational requirements of our PIC simulations 
the propagation of the pump pulse will not exceed 1 cm. 
We tend to avoid self-injection into the wake by adjust- 
ing the pump pulse intensity and the electronic density. 
Many different combinations of polarizations can be cho- 
sen for both waves, each of these possibilities results in 
particular force acting on the plasma electrons, when the 
two waves collide ^o| . The first point of this article is to 
summarize the dependance to pump and colliding pulse 
intensities, and to plasma density of this force. The rel- 
ative influence of stochastic heating and beat wave force 
on the injection mechanism, and later on the beam qual- 
ity will be discussed. After the choice of polarization in 
the blowout regime is clarified, we focus on the study of 
wakefield acceleration in the presence of an external, ho- 
mogenous, magnetic field and study its influence through 
simulations. The mechanisms leading to the enhance- 
ment of the beam quality will be examined. In a third 
part, fine-tuning of the counter-propagating low intensity 
pulse will be considered in order to limit the intensity of 
the external field. This situation will be illustrated in 
the case when the intensity of the pump pulse is raised 
to ao = 10. Then we will conclude. 
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II. SENSITIVITY OF BEAM INJECTION TO 
WAVE INTENSITY, PLASMA DENSITY AND 
POLARISATION OF LASERS 

A. Basic principles 

The Wakefield propagates in the plasma at the group 
velocity of the laser /3g defined by : /3g = Vg/c = 

( 1 — ^ j , where c is the speed of light, cjp and cjq re- 
spectively denote the plasma and laser frequencies. We 
use the quasi-static approximation and assume that the 
potential ^ created by the pump pulse only depends on 
= x—/3gt, where x and t denote the space and time coor- 
dinates normalized by c/cuo and uJq^ respectively. Then, 
the hamiltonian of an electron in the wakefield potential 
0, created by the pump pulse, reads : 



7l 



PgPx - HO, 



(1) 



where 7^ = 1 + p^. The normalized transverse mo- 
mentum is defined by p± = where is the trans- 
verse momentum and mg denotes the electron mass. The 
hamiltonian of a particle is an invariant H{^^px) = i^o, 
using ^ we deduce two solutions for the longitudinal 
momentum : 

Pt iO = MliHo + ± 7. {ll{H, + <t>f - ^If'^ , (2) 

where 7^ = 1 -\- /3g. The separatrix between trapped 
and untrapped orbits is given by a critical value, Hs = 
— — ^min, of the Hamiltonian. Replacing Hq by Hs in 
([2]) and retaining we get the two branches of 

Let us now comment the dynamics of an electron 
in vacuum, in the presence of two counter-propagating 
laser pulses. The hamiltonian reads H{ux^x) = 



(1 + . 



.1/2 



7, the longitudinal force acting on 



the electron is given by 

dUx dH 1 du^ 
dt dx 27 dx 



(3) 



Taking only into account the infiuence of the lasers, de- 
noted by their potential vectors Aq and Ai, we have 
= Ao + Ai. We will consider parallel linear polariza- 
tion (P linear polarization) and positive circular polariza- 
tion, that is u± = aocos{uJot — kox)ey-\-aicos{uJot-\-kox)ey 

and = ^cos(a;o^ — ^0^) + ^cos(cjo^ + ^0^) 

■^sm{ujot — kox) + ^sin(cJo^ + ^0^) respectively. 
Substituting the above expressions in Eq. (|3|) yields 

Fx = [alkosm{2{uJot - k^x)) 
Z7 

—a\kQ^\ii{2{ujQt + /cqx))] + ^o^o^i gi]2(2/cQx), (4) 



for the P linear polarization case, and 

Fx sm(2A:ox) ^5^, 
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(5) 



for the positive circular polarization case. When P linear 
or positive circular polarizations are used Eqs. (|4]) and 
(|5j) show the existence of a force F^^ spatially oscillating 
with a Ao/2 period. This force is not time-dependant, 
it is usually interpreted as ^ponderomotive force asso- 
ciated with the beat- wave jloMlSj. We should compare 
the beatwave force F^^^ to the longitudinal ponderomo- 
tive force, this latter scales as Fp = {2'y)~^ aQ/{cr) where 
r is the pulse duration. Taking the maximum value of 
Fi)u} the ratio Fi^y^/Fp becomes 



Fbw/Fp = 2aiUJor/ao. 



(6) 



When Fiju)/Fp > 1 electrons are trapped inside the Ao/2- 
long beatwave buckets and cannot be wiped out by the 
longitudinal ponderomotive force. On the whole for both 
polarizations electrons will undergo the beatwave force, 
therefore when the separatrix is such that mm{ux) < 
electrons will be trapped in the wakefield as a bunch Q • 
However in the P linear polarization case other terms are 
added to the force and the equations of the motion are 
no longer integrable, electron trajectories become chaotic 
[155. This phenomenon known as stochastic heating can 
provide very large momenta to some electrons [16-22|. 
When mm{ux) > the beatwave force is not efficient 
and electrons can hardly be trapped with positive circu- 
lar polarizations. In this case P linearly polarized collid- 
ing waves are necessary to give electrons the appropriate 
momentum in order to fill the gap between trapped and 
untrapped orbits. In the following two subsections, the 
sensibility of the injection process to polarization and in- 
tensity of the waves will be summarized. 



B. Low density: Strong dependance to polarization 
and wave intensity 

When the plasma density is about ~ 10~^nc a weak 
variation of the pump pulse is sufficient to modify the 
mechanisms allowing trapping of a bucket in the wake- 
field. To illustrate this strong dependance to laser pa- 
rameters we launched three sets of simulations with the 
PIC code C ALDER The simulation setup consists 
in two 30 fs linearly polarized waves with wavelength 
A = 0.8/im having their electric fields either linearly or 
circularly polarized (orthogonal linear polarization is de- 
noted by S). The pump pulse, which creates the acceler- 
ating wakefield, is focused to an 18/im full width at half 
maximum (fwhm). The peak normalized vector potential 
of the main pulse is associated to the laser intensity / by 
the formula a = 0.853 x 10-^A(/im)[/(W cm-^)]i/2, we 
considered a = 1.5 and a = 2. The low intensity pulse 
is counter-propagating and is focused to a 31 jam focal 
spot at a peak normalized vector potential ai =0.1 or 
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0.4. The waves interact with a mm-size plasma with a 
density = 4.3 x 10~^nc. The fluctuations between the 
two regimes are illustrated by Figs. [H When a = 1.5 and 
ai = 0.4 the ID separatrix between trapped and opened 
orbits is higher than = 0, in this case the beatwave 
force cannot provide enough momentum to electrons to 
push them in the wakefield. Nevertheless stochastic heat- 
ing due to P linear polarizations is a way to bridge the gap 
and inject a bunch in the wakefield (Fig. l(a)[ ). When 
a = 2 the ID separatrix is lowered (min(p^) < 0), in 
this case the beatwave force is enough to trap electrons 
in the wake therefore we can accelerate a beam using ei- 
ther P linear or positive circular polarizations. Note that 
no trapping occurs with negative circular polarizations, 
which is consistent with the theory introduced in subsec- 
tion III A[ Indeed straightforward algebra gives Fx = 0^ 
then the ponderomotive force hinders trapping of elec- 
trons. The quality of the accelerated beam will depend on 
the force which dominates during collision of the pulses. 
A relatively high intensity of the counter propagating 
laser (ai = 0.4) will foster stochastic heating i22|] 
and a higher charge will be injected with P linear polar- 
izations compared to positive circular polarizations (Fig. 
1(b)). On the contrary, if we reduce the intensity of the 
counter propagating pulse (ai =0.1) the beatwave force 



is favored and rules the injection mechanism (Fig. l(c)[ ). 
In this latter case, the beam quality is better without 
stochastic heating. 



C. Very low density: Weak dependance to laser 
polarization 

We now chose parameters relevant for the study of the 
blowout regime [gI, l23 - [26j . The plasma density was set 
to He = 2.5 X 10~^nc and the pump pulse intensity to 
a = the rest of the simulation setup was not modi- 
fied. Considering two circularly polarized waves rotating 
in the same direction (positive circular polarizations) is 
as efficient as considering P linear polarizations (Fig. [2]). 
It means that in this case cold injection [14, 27] due to 
the beating force is the key mechanism governing elec- 
tron injection. Stochastic acceleration weakly changes 
the number and the energy of electrons trapped in the 
wake field. The electron momentum distribution along 
the direction of propagation of the waves during their 
collision shows that the effects on electron dynamics of 
the two polarizations are close, as a result the electron 
energy distributions are almost identical (Fig. [2]). This 
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FIG. 1: Electron energy distribution obtained long after the 
collision of the two waves with rie — 4.3 x 10~^ nc. P lin- 
ear polarizations (blue curve). Positive circular polarizations 
(green curve). S linear polarizations (red curve) and negative 
circular polarizations (black curve, overlaid with axis in this 
case) . 



FIG. 2: Electron energy distribution from 2D PIC simula- 
tions, a = 4, ai = 0.1 and rie — 2.5 x 10~^nc. (a): Positive 
circular polarizations, (b): P linear polarizations, (c): Nega- 
tive circular polarizations, (d): S linear polarizations. 

dominant inffuence of the beating force upon injection 
was not clearly stated by Davoine et al. [14], here we 
underline that no relevant difference is triggered through 
use of P linear polarizations or positive circular polariza- 
tions in the blowout regime. The next section is devoted 
to the study of a new means to enhance beam quality 
after colliding pulse injection, in the blowout regime. 



III. INFLUENCE OF A MAGNETIC GUIDE 
FIELD 

The inffuence of a constant homogeneous guide mag- 
netic ffeld on LWFA is studied in this part. This mag- 
netic ffeld is assumed to be parallel to the direction of 
propagation of the waves. Electrons are still externally 
injected using a colliding counter propagating laser pulse 
[13]. The idea is to guide the electrons in order to im- 
prove the quality of the beam which is trapped in the 
wake ffeld. This mechanism was ffrst proposed in the 
context of LWFA, with a single pump laser (a = 3.5) 
and an electronic density rie = 3 x 10~^ nc prone to self- 
injection into the wake ffeld [28], with these parameters 



4 



self-injection can be dramatically enhanced and beam 
quality degraded. Here we aim at studying the influ- 
ence of a magnetic guide field in the blowout regime [sl 
with colliding pulse injection of the electrons, hence we 
chose a = 4 and Ue = 4.4 x lO^'^cm"^ = 2.5 x 10~^nc, 
thus abiding by ao > 4 and 2 < ao < 2uJo/ujp criteria 
proposed by Martins et al. p5|. The required magnetic 
field necessary to curve electron trajectories is about a 
hundred teslas [28j, such values are particularly strong 
but still available from the current pulsed magnet tech- 
nology [29I, the most advanced magnets can reach 90T 
for tens of ms durations and centimeter size lengths [30|] . 
The simulation setup consists in two 30 fs linearly po- 
larized counterpropagating waves with A = 0.8/im wave- 
length. They propagate along a constant homogeneous 
guide field Bq in a cm-long plasma, the normalised value 
of Bq is given by Bq — eB^/ {meOOo). Their electric fields 
are in the same plane (P linear polarizations). The pump 
pulse, which creates the accelerating wakefield, is focused 
to an 18/im full width at half maximum. The peak nor- 
malized vector potential for this pulse is still a = 4. The 
low intensity pulse is focused to a 31/im focal spot at a 
peak normalized vector potential ai = 0.1. The plasma 
frequency in the presence of a magnetic field can be ap- 
proximated by cj^ = (cjg+r^^)^/^ where uJrn andcjg repre- 
sent the frequencies of the magnetized and unmagnetized 
plasma, respectively. The cyclotron frequency is defined 
hy n = eBo/m. When Bq = 125T and Bq = 250T, one 
has l^^/cjg = 0.34 and l^^/cjg = 1.38 respectively which 
proves that the magnetization of the plasma may have 
some effect on the deformation of the wakefield as will 
be further discussed. Before the collision of the waves 
the magnetization has no infiuence on self-injection, no 
electron is trapped in the wake. 



Electronic density and transverse currents 
induced at the rim of the bubble 



Let us first identify the differences brought by the ad- 
dition of a magnetic guide field to the electronic distri- 
bution at the vicinity of the bubble boundaries (Figs. 
[3]). After the electron beam injection into the wake field, 
the strong ponderomotive force due to the main pulse 
still repels electrons and thus provides them longitudi- 
nal and transverse momenta. In the absence of magnetic 
field, electrons are submitted to the recall electric field 
induced by the bubble but the balance between this force 
and the ponderomotive force is favorable to the latter ; 
as a result electrons fiee along straight line trajectories 
(Fig. 3(a) ). When a longitudinal magnetic field is added 
electrons start to revolve around the bubble as a result of 
the magnetic force, this force combined with the electric 
force induced by the bubble completely modifies the dy- 
namics of the electrons. The gyro-radius of the electrons 
with py ^ is reduced when Bq is raised. Therefore the 
corresponding fiight path in the (x, y) plane appear to be 
bent, the trajectory will be even more curved when the 
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FIG. 3: Electron density long after the collision of the two 
waves. P linear polarizations. a = 4, ai = 0.1 and Ue = 



applied field is stronger (Figs. 3(b)|3(c)] ). When elec- 
trons revolve around the bubble they create a current 
perpend icularly to the plane of the figure, as evidenced 
by Figs. 4(a)||4(b)] Given the value of the electronic den- 
sity at the rear of the bubble (lxlO~^ < rie < 1.5x10"^), 
we deduce from Figs. [H the time averaged speed vf^ of 
electrons revolving around the bottleneck of the bubble 
0.5 < vf^^/c < 0.8. This current will act as a small 
solenoid, and thus the longitudinal magnetic field will be 
amplified. This feature and its influence on the beam 
dynamics will be detailled in the next subsection. 
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FIG. 4: Transverse component of the current density (nor- 
malized by ericc), the wakefield propagates in a magnetized 
plasma (B — 250 T). a = 4, ai =0.1 and rie 



B. A new mechanism to enhance the beam quality 

The intensity of the magnetic field is almost doubled 
locally (Figs. [5j) compared to the initial {t = 0) uniform 
map of Bx. We shall underline that this pattern is stable 
as we obtain quasi identical maps of B^ in this region 
of the bubble when the pump pulse has just entered the 
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plasma around ujot = 2280. Moreover we note that the 
geometry of the magnetic field lines is weakly altered by 
the electronic density modulations induced by the prop- 
agating bubble. Magnetic field lines stay almost parallel 
to the propagation direction. Let us now examine the 
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FIG. 5: Longitudinal component of the magnetic field (nor- 
malized by rUeOJo / e) , the wakefield propagates in a magnetized 
plasma Bq ^ 250 T {i.e. Bq 18.7 x 10"^). a = 4, ai = 0.1 
and ne = 2.5xl0~^nc. Magnetic field lines are superimposed 
(red curves) on Ex colormap. 

effect of the magnetic field on the dynamics of the ac- 
celerated beam. To evidence the differences between the 
magnetized and the unmagnetized regimes, we plotted 
kinetic energy density maps showing the evolution of the 
trapped beam at the rear of the bubble. In the unmag- 
netized case (Figs. [6]), the beam alternatively explodes 
(due to space charge effects) and focalises (due to the 
focalizing effect of the transverse electric field). This be- 
havior has a typical ~ 2000cj^^ period. Note that the 
beam acceleration is degraded because the components 
of the splitted beam will not see the maximal value of 
the longitudinal electric field (which is located on axis). 
The dynamics is completely different in the magnetized 
case (Figs. [7j), the longitudinal magnetic field is strong 
enough to curve the trajectories and hinder the explosion 
of the beam. As a result the beam is almost concentrated 
on axis, the transverse emittance is reduced and the main 
part of the beam, which also corresponds to the region 
where the magnetic field is the strongest, always sees 
the maximum value of the electric field. Next we will 
quantify the enhancement of the beam quality through 
the evolution of the energy distribution functions of the 
beam. 







- 50 
■ 0.3 

lo.25 






0.2 


1 




0.15 1 ° 
0., 






■o.os 


r r 


7400 7420 7440 7460 7480 


7700 7720 7740 7760 7780 



(c) 



(d) 



FIG. 6: Electron kinetic energy density (normalized by 
meC^ric), the wakefield propagates in a non magnetized 
plasma (Bq = 0). a = 4, ai = 0.1 and Ue — 2.5 x lO"'^ ric. 
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FIG. 7: Electron kinetic energy density (normalized by 
meC^ric), the wakefield propagates in a magnetized plasma 
{Bo = 250 T). a = 4, ai = 0.1 and Ue = 2.5 x 10"^nc. 



C. Enhancement of the beam quality 

As already mentionned, in the magnetized case, the 
beam is submitted to a more uniform accelerating field, 
this pattern boosts the particle acceleration leading to a 
slightly higher maximum kinetic energy when compared 
to the unmagnetized case (Fig. [8]). 

The focalizing magnetic field reduces the low energy 
tail of the energy spectrum, as can be seen in Figs 
[HI Obviously, this trend is enhanced when the guide 
field rises. With no guide field, the relative variation 



of the energy at full width at half maximum (fwhm) 
AEfwhrn/ Emax ^ 1% is excellent, but the rms value of 
the energy spread has small variations and reaches 7% 
at the end of the simulation (Table H]). When Bq = 125 
T, on the one hand the spread of the low energy tail of 
the distribution is reduced as shown by Fig. 8(a)[ and 
confirmed by the rms value ~ 4%, but on the other hand 
^Ef^hm/Emax IS slightly degraded. When Bq = 250 T, 
untrapped electrons carrying energies about 10 Mev con- 
centrate (ng locally reaches 1.5 x 10~^nc) at the rear of 
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uot 7000 9000 11000 13000 

^0 (T) 

5.33% 5.88% 4.53% 6.61% 

(2.05%) (1.09%) (0.81%) (0.97%) 

125 5.50% 5.45% 3.41% 3.96% 

(5.03%) (4.89%) (1.94%) (1.58%) 

250 3.96% 2.67% 2.84% 2.91% 

(3.72%) (2.75%) (1.54%) (0.63%) 



TABLE I: Evolution of the half width of the electron dis- 
tribution function with a = 4, ai = 0.1, rie = 2.5 x lO"'^. 
Standard deviation, i.e. rms value, and relative variation 
{AEfwhm/Emax) of the distribution, where the subscript 
fwhm denotes the Full Width at Half Maximum. 




(a)Bo = 125 T (bold lines). (h)Bo = 250 T (bold lines). 



longitudinal gradient of the electronic density to enhance 
trapping [1^,^,[32|. In an alternate approach, assuming 
initially homogenous plasma, one can slowly evolve the 
laser pulse shape to alternate periods of expansion and 
contraction of the bubble, to respectively trigger and stop 
self-injection [33| of the electrons into the bubble. How- 
ever this technique seems hard to adjust to get a unique 
mono-energetic bunch. In this paper, injected electrons 
are confined by using a magnetic guide field but we note 
that the intensity of the field, required to substantially 
enhance the beam quality, depends on lasers pulse shapes 
and durations. In this section we did not pay attention 
to the tuning of the low intensity colliding pulse, we have 
shown that the guiding induced by Bq is enough beyond 
some threshold, function of the parameters of the simula- 
tion. However, if we decide to lower the blowout stability 
by increasing the main pulse intensity, the required in- 
tensity of the guide field grows to values largely out of 
reach of the current technology, then fine-tuning of the 
colliding pulse becomes necessary. The next section is 
devoted to this issue. 



IV. INFLUENCE OF A MAGNETIC GUIDE 
FIELD AT HIGHER INTENSITIES 



FIG. 8: Electron energy distribution from 2D PIC simu- 
lations. P linear polarizations. a = 4, ai = 0.1 and 



2.5 X 10" 



Dashed lines correspond to B = 0. 



the bubble. These low energy (i.e. < Ek < 25 Mev) 
electrons are evidenced by bumps in the beam energy dis- 
tribution (Fig |8(b)] ). This low energy bump slowly slides 
out of the simulation box as these electrons are not in- 
jected in the wakefield, and therefore should not be con- 
sidered for the interpretation of the diagnostics concern- 
ing the accelerated beam. According to this comment, we 
note that a 250 T guide field is enough to completely sup- 
press the low energy tail during the whole simulation. A 
clear enhancement of the beam quality is obtained, first 
the final rms value is below 3% and AEf^hrn/ Emax < 1% 
(table H]) thus providing a very sharp control on the final 
energy of the beam, and second the number of electrons 
at the highest energies does not vanish, as in the unmag- 
netized case, but on the contrary grows up to 60 x 10^ 
part/Mev, nearly twice the value of the unmagnetized 
case ! To our knowledge such an acute mono-energetic 
electron beam production, with complete extinction of 
the low energy tail has never been evidenced. 

For a given pump pulse, two well-known controllers of 
the low intensity laser can be used to optimize the beam 
quality. The low intensity pulse duration can be mon- 
itored and the transverse fwhm of the spatial envelope 
adjusted, these two parameters together usually make it 
possible to get a quasi-mono-energetic electron beam in 
the blow-out regime [1^ l25| . There are other ways to en- 
hance the beam quality. For example, one can resort to a 



A very high intensity wave is considered now, the pump 
pulse which is linearly polarized is assumed to have a 
peak normalized intensity a = 10, and a duration of 30 
fs. It is focused to a 35.65/im large focal spot. The col- 
liding pulse has a peak normalized intensity ai = 0.1, 
a 30fs duration and three focal spot sizes were consid- 
ered D = 60/im, 36 /im and 10 /im. The wavelength 
of the waves is A = 0.8/im. In these simulations, the 



plasma density is still rie = 2.5 x 10 



Unless other- 



wise mentionned the following simulations were run with 
Bq = 125T. When considering a = 10 and a counter- 
propagating wave focused to a 10 /im focal spot, one has 
a paramount effect of the magnetic field on the distri- 
bution function (Table HIl Fig. [9]). The electron energy 
distribution becomes almost mono-energetic. After the 
laser pulse has propagated through the plasma by 3.8 
mm, the electron energy distribution is still quite mono- 
kinetic and the maximum electronic energy exceeds one 



GeV (Fig. |9(b) ). Table HIl shows the evolution of the en- 
ergy spread of the electron energy distribution with time. 
When no field is applied the quality decreases whereas 
we get an acute control on the beam energy with Bq. It 
must be pointed out that the accelerated charge is small 
(close to 50 picocoulombs). The electron energy distri- 
butions corresponding to four focal spots are compared 
at some time. Figure [13] shows that the distribution be- 
comes much more mono-energetic when the focal spot of 
the perturbing wave is smaller than the one of the main 
pulse. The magnetic field is much more efficient when 
considering a small value of D. Then, we have checked 
that this peak exists in a very small range of values close 
to ai = 0.1. When ai = 0.08 the distribution function 
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TABLE II: Evolution of the half width of the electron dis- 
tribution function, a = 10, ai = 0.1, Ue = 2.5 x 10~^. 
Standard deviation, i.e. rms value, and relative variation 
(AEfujhm/Emax) of the distribution. 



shows a lower magnitude peak, and the charge acceler- 
ated in the first bubble is about 10 picocoulombs. The 
peak still exists when ai = 0.102. But in the case of 
a higher value of ai(for instance ai = 0.15) the energy 
distribution does not show a very thin high-energy peak 
any longer because the injected charge is higher and the 
magnetic field is too weak to concentrate the beam effi- 
ciently. 
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nitude of the magnetic guide field is increased to 250 T, 
accordingly with the results obtained with a lower inten- 
sity of the pump pulse (ao = 4) in section ITlI CI Figure [12] 
shows that, a shorter pulse duration {At = 10 fs) for the 
counterpropagating wave, also makes the electron energy 
distribution more mono-energetic. Then, in order to ob- 
tain a more mono-energetic distribution when a = 10 and 
ai = 1, a very strong magnetic field Bq = 250T and a 
short duration for the counterpropagating wave At = 10 
fs were considered (Fig. [T3|) . As expected, a high en- 
ergy peak is obtained. One should point out that the 
charge accelerated in the first bubble is still close to 50 
picocoulombs. To summarize when ai is too high, that is 
when the charge injected in the wakefield exceeds some 
treshold, we can not stop a drop of the beam quality by 
imposing an external field alone, at least the main pa- 
rameters (focal spot size and duration) of the colliding 
pulse shall be reduced. 




(b)cjo^ = 30000 
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FIG. 11: Electron energy distribution from 2D PIC simula- 
tions at (jJot — 13000 with a — 10, ai — 1^ rie — 2.5 x 10~^nc, 
D = 10/im and At = 30fs. (a)5o = 0, (b)Bo = 125 T, 
(c)Bo = 250T. 



FIG. 9: Electron energy distribution from 2D PIC simulations 
with a — 10, ai — 0.1, rie — 2.5 x 10~^nc and D — lO/xm. 
The dashed line corresponds to Bo = and the bold one to 
Bq 125 T. 
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FIG. 10: Electron energy distribution from 2D PIC simula- 
tions at cjot = 13000 with a = 10, ai = 0.1, rie = 2.5 x 10~^nc 
and Bq = 125 T. {s.)D = 59.5/im, {h)D = 35.7/im, {c)D = 
20/im, (d)^ 10/im. 

No high-energy peak was seen in the distribution func- 
tion for larger values of ai, actually the same kind of 
distribution is obtained when ai = 0.2, ai = 0.5 and 
ai = 1. Figure [TTI shows that the electron energy dis- 
tribution becomes more mono-energetic when the mag- 




Ek (MeV) 

FIG. 12: Electron energy distribution from 2D PIC simula- 
tions at cjot — 13000 with a — 10, ai — 1^ Ue — 2.5 x 10~^nc, 
Bq = 125 T and D = lO/xm. (a) At = 30fs, (b)At = lOfs. 



V. CONCLUSIONS 

In the first part of this paper, we summarized some 
essential results about the most efficient choice of po- 
larizations for injection in the bubble, in the colliding 
pulse scheme. At rather high electron density (^ 10~^nc) 
and moderately relativistic electromagnetic wave inten- 
sity (ao ^ 2), more particles are accelerated to high 
energies in the case of P linear polarizations that is 
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E^MeV) 

FIG. 13: Electron energy distribution from 2D PIC simula- 
tions at LOot — 30000 with a — 10, ai — 1^ rie — 2.5 x 10~^nc 
and Bq — 250 T. The spatio-temporal shape of the laser en- 
veloppe is defined by D = lO/xm and At = lOfs. 

to say when electrons undergo the action of the beat- 
wave force and all the others. For higher intensities and 
lower densities (~ 10~^nc) the beatwave force, that is 
cold injection, can be more efficient. The second and 
main part of this paper has been devoted to the study 
of the influence of an external static magnetic field on 
the Wakefield acceleration process, within the colliding 
pulse scheme. To our knowledge this idea was never ex- 
plored. The magnetic field is supposed parallel to the 



direction of propagation of the two counter-propagating 
waves. It has been shown that the Bq field creates a 
transverse current, the latter current can induce a raise 
of Bx at the rear bottleneck of the bubble. Therefore 
the beam dynamics is substantially modified as the beam 
is constrained to stay in the maximum acceleration re- 
gion of the bubble. Beam emittance is considerably 
reduced and maximum kinetic energy slightly boosted 
compared to the unmagnetized case. This mechanism 
provides means to dramatically enhance the beam qual- 
ity in the blowout regime. We achieved tremendous 
amelioration with the setup: a = 10, ai = 0.1 and 
He = 2.5 X 10~^nc. After roughly 4 mm of wakefield 
acceleration, without Bq field the electronic energy dis- 
tribution is noisy AEf^hrn/ Emax ^ 100% whereas we 
get AEfwhm/Emax ^ 3% whcu the plasma is magne- 
tized with a 125 T field. Nevertheless the intensity of the 
^0 field may be limited by technological considerations 
[30], thus acute control of the beam quality may require 
some fine-tuning of the colliding pulse parameters. For a 
given pump pulse, one should adapt the intensity, dura- 
tion and focal spot size of the counter-propagating laser 
pulse. 
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